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ABSTRACT 

We examine stellar population gradients in ^ 100 massive early type gala xies spanning 1 80 < < 370 km s“* 

and Mfc of -22.5 to -26.5 mag, observed as part of the MASSIVE survey (iMa et alj|2014l) . Using integral-field 
spectroscopy from the Mitchell Spectrograph on the 2.7m telescope at McDonald Observatory, we create stacked 
spectra as a function of radius for galaxies binned by their stellar velocity dispersion, stellar mass, and group 
richness. With excellent sampling at the highest stellar mass, we examine radial trends in stellar population 
properties extending to beyond twice the effective radius (~ 2.5Re). Specifically, we examine trends in age, 
metallicity, and abundance ratios of Mg, C, N, and Ca, and discuss the implications for star formation histories 
and elemental yields. At a fixed physical radius of 3-6 kpc (the likely size of the galaxy cores formed at high 
redshift) stellar age and [a/Fe] increase with increasing cr* and depend only weakly on stellar mass, as we might 
expect if denser galaxies form their central cores earlier and faster. If we instead focus on 1 - 1 .5Re, the trends in 
abundance and abundance ratio are washed out, as might be expected if the stars at large radius were accreted by 
smaller galaxies. Finally, we show that when controlling for cr*, there are only very subtle differences in stellar 
population properties or gradients as a function of group richness; even at large radius internal properties matter 
more than environment in determining star formation history. 

Subject headings: galaxies: elliptical and lenticular, cD, galaxies: evolution, galaxies: kinematics and dynamics, 
galaxies: stellar content 


1. INTRODUCTION 

The assembly history of elliptical galaxies remains a major 
unsolved problem for galaxy evolution. Recent observations 
point to dramatic size evolution of the most massive galax¬ 
ies from z ft: 2 to the p r esent (e.g., lyanDpkkumeL^ 200^ 
van der Wei et al.l 120081: iPatel et 5l 2013 : Ivan der Wei et al.l 


spite more than thirty years of effort, most observations of 


2014 . The extent t o which these tre nds require late-stage 


minor mergers (e.g., lOser et afl l2012h or can be explained 
by the addition o f large r, younger galaxies a t later times 
( Valentinuzzi et al.1 1201 Ot iNewmian et akl 1201 2t iBarro et aH 
1201 3h remains a topic of ongoing debate. Information lurking in 
the faint outer parts of present-day massiye ellipticals can com¬ 
plement high-redshift measurements. Radial gradients in stellar 
populations distin guish when and how the stars at large radius 
were f ormed (e.g., Whitd 1980HKobayashill2004lGreene et al.1 
1201 31: iHirschmann et al.M2015li . while the kinematics of the 
stars (e.g., V/ct*, the leyel of radial anisotropy, etc. ) contain 
clues about how these stars entered the halo (e.g., Wu et al.l 


i 


20l4lArnoldetal.l20l4lRaskuttietal.l20l4lNaabetal.l201 

Here we focus on the ayerage radial trends in stellar pop¬ 
ulations of early-type galaxies using our ambitious suryey 
of the hu ndred most MA SSIVE galaxies within one hun¬ 
dred Mpc (iMa et al.ll20l4 . Obseryations of the stellar pop¬ 
ulations in elliptical galaxy outskirts are challenging, since 
their surface brightnesses drop steeply with radius. De- 


half-light radius (15 

spinrad & Tayloi 19711: Faber etal. 1977: 

Gorgasetal.l 199C 

Fisher etal.1 19951: Kobayashi & Arimoto 

1999: 

Ogando et al. 20051: Brough et al.1 20071: iBaes et al. 

2007, 

Annibali et al.l 20071: Sanchez-Blazo 

uez et al.l 20071: 

Rawle et al. 20081: Kuntschner et al.l 20101 

McDermid et al. 

20151: Oliya-Altamirano et al.ll2015l). Resolyed stellar nonula- 


tion studies haye uncoyered a low-metallicity halo component 
at yery l arge radius, but only in a handful of nearby galax¬ 
ies (e.g., Kaliraietal. 2006t Harris et alJ 119991: iReikubaet al 


2005; lHams_et^ 20071 ICrnoieyic et~ 2013HPastorello etaT.. 

2014 : Peacock et al.1 2^51; I Williams et al. 2015 ). There are 


(Carollo et al. 1993 

: Carollo & Danziger 1994 

: Mehlert et al.1 

20031 Kelson et al.l 

20061 LSnolaor etal.1120101: 

Pu et al .1120lot 

Pu&Han 20111). 

Eyen integral-field spectrographs, now 


widel y used for the study of spatially resolyed galaxy prop¬ 
erties ( Emsellem et al.ll2004 ISarzi et a^l.ll2005 ICappellari et al.1 


l2006l 2012 ). include few obseryation s that extend beyond 
the half-light radiu s in integrated light (IWeiimans et al.l 120091: 
iMumhy et al.ll20Tlh . 

The ^ 100 galaxies analyzed here represent a significant im- 
proyement oyer preyious work. Using coadded spectra as a 
function of radius, we will inyestigate whether radial gradients 
depend not only on stellar yelocity dispersion, but also on other 
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intrinsic galaxy properties such as stellar mass and environmen¬ 
tal density. We find a hint that at fixed cr*, galaxies residing in 
higher densities are older and more a-element enhanced, but 
we see no evidence for differing gradients for galaxies in low 
and high densities. 

We present the sample in ^ the observations and data reduc¬ 
tion in ^ our stellar population analysis in ® and the radial 
variations in stellar populations in ^ We summarize our find¬ 
ings and conclude in ^6] Throughout we assume a concordance 
cosmology with //p = 70 km s“^ Mpc“', Hm = 0.3 and = 0.7 
(iDunklev et al.l[200^ . 


2. SAMPLE 

The MASSIV E sample selection is described in detail in 
iMa et al.l (|2014|) . Eor completeness, we summarize the sam¬ 
ple selection briefly here. MASSIVE is a volume-limited sam¬ 
ple of the 116 most massive galaxies within 108 Mpc. The 
galaxies are selecte d from the 2MASS Redshift Survey (2MRS; 
iHuchra et al.ll20l3) , using a total /T-band magnitude limit of 
Mk < -2 5.3 mag (roughly 1 0*^ ^ Mq). Using the Hyperleda 
database (iPaturel et alJ 1200^ . we apply a morphological cut 
to remove large spiral and interacting galaxies. The result¬ 
ing 116 galaxies span a wide range of stellar velocity disper- 
sion 180-400 km s“') and, based on the group catalog of 
ICrooketalJ(l2007h . are found in a wide array of environments, 
from “field” galaxies with no ^ L* companions to rich clusters 
(Coma, Perseus, and Virgo). In this paper, we include 49 MAS¬ 
SIVE galaxies with large-format integral-field spectrograph ob¬ 
servations in hand. While we are also obtaining high spatial res¬ 
olution integral field spectroscopy for a subset of galaxies fo r 
black hole mass determinations (e.g.. McConnell et al.11^1^ . 
this paper deals exclusively with the Mitchell data. 

In addition to the ongoing MAS SIVE survey, we inc lude here 
a sample of lower-mass galaxies (iGreene et alj|201^ selected 
direct ly from the Sloan Digital Sky Survey (SDSS; York et all 
l2000h and observed with an identical set-up with the Mitchell 
Spectrograph. As th e selection was different, we briefly re - 
view it here (see also lGreene et al.lF2012l : lRaskutti et al.ll2014l) . 
Since the spectral resolution of the Mitchell Spectrograph is 
(Tinst ~ 150 km s“' at 4000A, we select galaxies with disper¬ 
sion measurements from the SDSS that are greater than this 
value. Individual fibers are A'.'2 in diameter, and so we aim 
for galaxies with effective radii at least twice as large. Galax¬ 
ies with distances of 40-95 Mpc are large enough to be well- 
resolved but small enough t o fit into one pointin g. We use a 
color selection of M-r > 2.2 dStrateva et al.ll200'll) . which pref¬ 
erentially selects early-type galaxies, and then remove the few 
edge-on disk galaxies by hand, but keep SO galaxies. These 
46 galaxies are of uniformly lower stellar mass than those in 
MASSIVE. Taken together, the sample galaxies span a range of 
-22.5 < Mk < -26.5 Vega mag. The sample properties for the 
95 galaxies (49 from MASSIVE) are summarized in Figure|2] 


3. OBSERVATIONS AND DATA REDUCTION 

The majority of the galaxies presented here were observed 
over ten observing runs between Sept 2011 and April 2014. 
There are an additional four gala xies taken earlier as part of the 
PhD thesis of Jeremy Murphy dMurphv et al.l 2Mlh a nd four 
taken by Nicholas McConnell (iMcConnell et alJl2oi2h . All of 
these data were taken with the George and Cynthia Mitchell 
Spectrograph (t he Mitchell Spectrograph, formerly VIRUS-P; 
iHill et alllioO^ on the 2.7m Harlan J. Smith telescope at Mc¬ 


Donald Observatory. The Mitchell Spectrograph is an integral- 
field spectrograph with 246 fibers subtending 472 each and cov¬ 
ering a 107" X107" field of view with a one-third filling factor. 
As such, the Mitchell Spectrograph is ideal for stud ying the low 
surface brightness outer pyts of nearby galaxie s (iBlanc et al 
2009; lYoachim et al.l 1201^ iMurphv et all 120111: lAdams et al 


2012 ). 


We utilize the blue setting of the Mitchell spectrograph with 
a resolution of R 850 and spanning a wavelength range 
of 3550-5850 A. This resolution (roughly a FWHM of 5A) 
delivers a dispersion of ^ 1.1 A pixel"' and corresponds to 
(T* Ri 150 km s"* at 4300 A, our bluest Lick index, improving to 
^100 km s"* at the red end of the spectrum. Each galaxy was 
observed for a total of ^ 2 hours on source with one-third of 
the time spent at each of three dither positions to fill the field of 
view. Observations are interleaved with sky observations of 10 
min duration. Initial data reduction is accomplished using th e 
custom code Vaccine (I Adams et al.ll201 ll; Murphy e t al.ll201 ll). 
The de tails of our data reduction are described in Murphy et al.l 
dloTih . so we repeat only a brief overview for completeness 
here. 

We first perform overscan and bias subtraction for all science 
and calibration frames. The fiber trace is determined from the 
twilight flats, taking into account cur vature in the sp atial di¬ 
rection and following the techniques of lKelsoiil (l2003l) to avoid 
interpolation and thus correlated errors. Arcs are used to de¬ 
rive a wavelength solution with typical rms residual variations 
about this best-fit fourth-order polynomial between 0.05 and 
0.1 A. The twilight flats are also used to construct the flat field, 
once the solar spectrum has been modeled and removed. The 
flat field is typically stable to <0.1 pixels for typical thermal 
variations in the instrument. The flat field is then applied to 
all of the science frames to correct variations in the pixel-to- 
pixel responses, as well as the relative fiber-to-fiber variation, 
and the cross-dispersion profile shape for every fiber. The sky 
is modeled using off-galaxy sky frames observed with a sky- 
object-object-sky pattern. The sky frames are processed in the 
same manner as the science frames. In general, each sky nod 
is weighted equally, except for very cloudy conditions. Finally, 
cosmic rays are identified and masked. 

We use software d eveloped for the VENGA project 
(iBlanc et al.ll2009l 1201 3h for flux calibration and final process¬ 
ing. We observe flux calibration stars using a six-point dither 
pattern and derive a relative flux calibration in the standard way. 
We can test the wavelength dependence of the flux calibration 
by comparing the shape of the spectrum in the central fiber of 
the Mitchell Spectrograph with the SDSS spectrum for those 
galaxies with SDSS spectra. We find < 10% disagreement in 
nearly all cases, with no more than ~ 15% differences at worst. 
We then correct the default astrometric solutions using pho¬ 
tometry, by de riving an astrome tric match betwe en each image 
and a n SDSS (lYork et al.ll2000l) or PanSTARR S dSchlaflv et all 
1201 2t iTonrv et al.ll2012l: Magnier et al.ll20r^ image. Finally, 
all fibers are interpolated onto the same wavelength scale and 
combined into radial bins. 


3.1. Effective radii 

In our previous work, we adopted the SDSS model radius 
(the de Vaucouleurs fit) as the effective radius (Re)- While 
there is considerable evidence that the shape of th e light pro¬ 
file changes systematically with galaxy mass (e.g ., ICaon et aU 
Il993t iFeirarese et al.ll200d iKorinendv et al.ll200^ . 









































































MASSIVE II 


3 


M* {M^) 

IQIO 1011 1018 




Fig. 1.— Distributions of basic structural properties for our sample. On bottom left, stellar velocity dispersion (cr*) measured within the central fiber is plotted 
against Mk (mag; total magnitude) from 2MASS. On top left, we show effective radius, from the NSA catalog (or from 2MASS with a correction; §3.2) against Mk 
and on bottom right effective radius versus stellar velocity dispersion. 



Fig. 2.— Stellar velocity dispersions (cr*) measured within the central 4'/2 
Mitchell fiber are compared with literature values. Open symbols have litera¬ 
ture measurements from the SDSS, while filled symbols are from Hyperleda. 
The overall agreement is quite good, with (cr*Ms—<7*LitV = 0-01 =b0.09. 

fitting the galaxies with a fixed Sersic index of four has the 
benefit that we are less sensitive to both sky su btraction errors 
dMandelbaum et aDl2005t iBernardi et alJ[200^ and to the de- 
tailed shape of the ligh t prohle in the very faint wings (e.g., 
iLackner & Gunnll2012l) . In the case of MASSIVE, - 1/3 of 
the galaxies do not have SDSS imaging, so we adapt a size mea¬ 
surement from 2MASS (iJarrett et al.ll2003h . We use the median 
effective radius measured from the JHK band (see Eq. 3 in 
iMa et alJiMl) . The 2MASS size measurement tends to un- 


deres t imate the galaxy sizes relative to the SDSS (iLauer et al.l 
l2007l iKormend v et al.ll2 009h. T o put the size measurements on 
equal footing, iMa et alJ (120141) fit a linear conversion between 
f^2MASS and RsdSS ■ logio(R2MASs) = 0.81ogjg(RsDSs) “ 0.076 
(their Eq. 4). We use this relation to correct the 2MASS mea¬ 
surements to match the SDSS measurements. 

3.2. Radial bins 

Spectra from individual hbers, with the exception of those 
at the very center of the lEU, have inadequate signal for stellar 
population studies. Therefore, all of our analysis is performed 
on binned spectra. We utilize two binning schemes here, in both 
cases dehning elliptical annuli based on the axis ratio measured 
by the SDSS or 2MASS. Since galaxies get larger as they get 
more massive, we make bins of width Q.SRg. However, it is 
also interesting to look at trends as a function of physical size, 
and here we extend from 0-15 kpc in 3 kpc increments. Eor 
the physical bins, we use the central fiber as the central bin; 
while this central fiber corresponds to a different physical size 
for each system, it provides our highest spatial resolution bin. 

3.3. Stellar velocity dispersion measurements 

Stellar velocity dispersions are required for measuring Lick 
indices. This is because at fixed intrinsic absorption, as the ve¬ 
locity dispersion increases, the measured EW decreases. Thus, 
a correction must be applied to put all indices on the same 
scale. We also use the central stellar velocity dispersions to 
rank galaxies in constructing coadded spectra. 
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Rest Wavelength (A) 

Fig. 3.— The coadded spectra within 0.5Re, for stellar velocity disper¬ 
sion stacks starting from cr* < 220 km s“^ (bottom) to cr* > 290 km (top). 
Each spectnim has been normalized to unity as described in the text, and the 
flux density offsets are arbitrarily applied for display purposes. In creating the 
stacks, each galaxy was smoothed to a value 30% more than the upper end of 
the dispersion range in creating the stack. 


Our central 4"2 fiber is similar to the SDSS fiber and thus eas¬ 
ily compared w ith the literature. We use pP XF to measure 
the dispersions dCappellari & Emsellemll2()0^ . We compare 
our (T* measurements with SDSS me asurements whe n avail¬ 
able or our compiled literatu re values ( Mae^^ l2014i) derived 


from the Hyperleda database dPaturel et al.ll2003h . We find very 


good agreement, with the median ((Tms “CrLitVcrLit = 0.01 ±0.09. 
That is, we see no systematic offset between the two sets of 
measurements and a scatter of 9%. We compare our mea- 
surements with the literat ure in Figure |2] As discussed in 
Ivan den Bosch et ^ d2015l) . there is some excess scatter at the 
high dispersion end, specifically when we compare with Hy¬ 
perleda. However, we hnd better agreement with Hyperleda 
than reported by van den Bosch et al. predominantly due to re¬ 
cent (post-2013) changes in Hyperleda. We have 51 objects 
in common with the van den Bosch HEX catalog, and we also 
hnd reasonable agreement with their cr* measurements, with 
(O'MS - o-het)/ o-het = -0.03 ± 0.07. Throughout, cr* refers to 
our measurements from the central hber unless otherwise spec- 
ihed. 


4. ANALYSIS 


4.1. Stellar population modeling approach 

We use Lick indices as a tool to trace the stellar popu¬ 
lations. Lick indices were developed as a way to extract 
stellar population information from spectra without Hux cali- 
brati on, but still circumvent classic age - metallicity degenera¬ 
cies dBursteinl 19^ iFaber et ^ 1 19851 : I Worthev et al.Tl 1 99^ 
iTrager et al.l 19981) . The Lick indices are narrow regions of the 
spectra (typically ^ 20A wide) that are dominated by a single 
element and thus are predominantly sensitive to orthogonal as¬ 
pects of the stellar population properties; H/3 is sensitive to age, 
Le indices to [Fe/H], and so on. It is important to bear in mind 
that at the velocity dispersion of our target galaxies (200-400 
km s“') all indices are blends of multiple ele ments. See, for 
example. Table 1 in iGraves & SchiavonI (l2008l) for the primary 


elements that dominate the Lick indices used in this work. 

Lick indices are still widely used in the literature as they 
mitigate difficultie s in modeling the ef f ects of abundance ra ¬ 
tio changes (e.g., Worthev et ^ 1 1994t iGallazzi et al.l l2005l) . 
iTrager et al.l (l2000blal) developed a technique to derive index 
responses from stellar atmosphere models, such that even if 
the full spectrum cannot be calculated, the index EW changes 
due to changing abundance ratios can be incorporated into an 
analysis of line EWs. On the other hand, the current genera- 
tion of full spectral synthesis codes are very sophisti cated (e.g., 
IVazdekis et al.l 1201 Ot IConrov & van DokkumI 1201 2|) and have 
been designed to fit non-solar abundance ratios. 

For reference, we review the main indices that we use to de¬ 
rive the basic stellar populat i on pa rameters. The reader is re¬ 
ferred to lGraves & Schiav onI (1200 8) for more de tailed informa¬ 
tion. We use lick ew JOraves & Schiavonll2008l) to measure the 
Lick indices and the stellar population modeling code EZ_Ages 
(iGraves & Schiavonll2008h to convert the Lick indices to phys¬ 
ical parameters (age, [Fe/H], [a/Fe]). The code works on a 
hierarchy of index pairs, starting with H/3 and (Fe), and itera- 



rectly measure the oxygen abundance and oxygen is the most 
abundant heavy element, it is misleading to quote total metal- 
licity. Instead, we quote [Fe/H], which is directly inferred 
from the Fe indices. If we assume that [O/Fe] tracks [Mg/Fe], 
then we can use the latter to infer [Z/H]. We will generally as¬ 
sume that O and Mg follow similar trends as they are both a 
elements, and thus use [Mg/Fe] interchangeably with [a/Fe]. 
Based on the same assumption , we also will use the conver¬ 
sion from ITrager et all ( l2000bl) : [Z/H]=[Fe/H]H-0.94[a/Fe] to 
calculate the metallicity. For alter nate approaches t o mod el¬ 
ing oxygen us i ng Lic k indices, see Ijohansson et al.l ( 1201 2|) or 
IWorthev et al.l (|2014|) . I n our default runs, w e utilize the a- 
enhanced isochrone from ISalasnich et ^(l2000t) and the default 
assumption that [0/Fe]= 0.5 to match the a-enhanced isochrone 
value. 

Carbon is roug hly half as abundant by number as oxygen at 
solar abundances (lAsplund et al.ll20 6^. We deriv e [C/Fe] from 
the C 7.4668 index (e.g Tripicco & Belli 19^ ITrager et'd] 
Il998t Graves et'aniioOTt iJohansson et alJl2012l) Because the 
amount of C locked into CO depends on the oxygen abun¬ 
dance, as we lower the assumed O abundance, the i nferred C 
abundance drops commensurately (IServen et al.l2005]) . We will 
quantify the magnitude of this effect below. 

Nitrogen (one-tenth the O abundance by number in the Sun) 
is then derived from CNl (41 43-4178A) once we have a C 
abundance. As discussed in iGreene et al.l (1201 3l) . there are 
some uncertainties associated with CN due to the low S/N at 
the blue end of the spectrum, but our nuclear CN measurements 
match those from the SDSS spectra of the same galaxies (our 
Mitchell indices are 0.01 ±0.02 mag lower than the SDSS mea¬ 
surements), giving us some conhdence in these measurements. 

Finally, the Ca abundance (one-hundredth the C abundance 
by number in the Sun) is based on the Ca4227 index (4223.5- 
4236.OA; IWorthevlll998h . This index is blended with CN, and 
thus the [Ca/Fe] measurement is most uncertain, as it is de¬ 
pendent on both the [C/Fe] and [N/Fe] measurements. Others 
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Fig. 4.— Radial gradients in age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] as calculated by EZ_Ages from the Lick indices measured in the coadded spectra. 
The measurements are made on four stacked spectra binned on stellar velocity dispersion (see figure key), and are shown as a function of R in kpc (left) or R/Re. 
We fit the radial gradients with a power law of the form X = A log(/?/R 3 ) + B for each stellar population parameter X, where R^ is either 3 — 6 kpc or 1 — \.5Re. The 
fits to the highest (long-dashed lines) and lowest dispersion (dotted lines) are shown here, and in Tables 1 & 2. Note the decline with radius in [Fe/H] and [C/Fe] in 
contrast with the radially constant age, [Mg/Fe], [N/Fe], and [Ca/Fe]. To indicate systematic errors in the light elements due to the unknown oxygen abundance, we 
also show the resulting models assuming [0/Fe]= 0.1 rather than the default [0/Fe]= 0.5 (keeping [0/Fe] constant with radius in both cases; a* > 290 km s~^ coadd 
in dark and 220 < cr* < 250 km s“^ in light color). The [C/Fe] lines with alternate oxygen abundance have been offset by -0.1 dex for presentation purposes. 
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have used Ca H+K (IServen et alJl2005l : IWorthev et alJl2014l) or 
the Ca Triplet at 8600A, vyhere the last is als o sensitive to the 
dwarf-to-giant ratio (e.g.. lCenarro et a02004l) . 

Since the Lick indices are very sensitive to small errors in sky 
subtraction and other small-scale errors, we construct stacked 
spectra and measure the average radial trends in the Lick in¬ 
dices and resulting stellar populations. 


4.2. Equivalent widths and emission line corrections 


There are a number of systematic effects that may impact the 
equivalent widths (EW). One relates to the velocity dispersion: 
while the default dispersion corrections within lick_ew were 
verified originally over the range cr* = 250-300 km s“', some 
of our galaxies exhibit even larger cr,. We therefore check our 
disper sion corrections using simple stellar population models 
from dConrov & van DokkumI 12013) . broadened over the full 
observed range of cr*. The Lick indices that we recover us¬ 
ing the default lick_ew corrections agree with the input values 
within 0.02A for an a-enhanced model with a Salpeter IME, 
which is a good approximation to our galaxies. The corrections 
are a weak function of stellar population parameters, but only 
at the hundredths of an A level, which is small compared to our 
other sources of systematic error. 

A larger correction must be made for low-level emission that 
can fill in the absorption lines and artificially lower their equiv¬ 
alent widths (EWs). Weak emission from warm io nized gas is 
very common in the center s of elliptical galaxies dSarzi et al.l 
l2010t lYan & Blantonll201^ . and small amounts of line infill 
can lead to significant errors in recovered parameters. Even 
0.1 A errors in Hff EW can lea d to errors of ^ 1-2 Gyr in the 
modeling (e.g.. lSchiavonl[200^ . 

Given the very low emission levels, and the large uncertain¬ 
ties involved, we compare two met hods for determining the 
levels of [O III] and H/3 emission. In iGreene et al.l ( 20121) we 
utilized pPXE-i-GANDAL F developed by M. Sarzi ( Sarzi et all 
|2006|) and M. Cappellari dCappellari & Emsellemll2004 ) to si¬ 
multaneously model the stellar absorption and emission lines. 
GANDALF is yery robust for well-detected line s, but under¬ 
constr ained for yery weak emission. Following iGreene et al.l 
d2013l) . we also fit each spectrum w ith an empirical tem plate 
drawn from the composite spectra of iGrayes et alJ d2010l) . We 
then fit the [O III] emission in the residual spectrum, and sub¬ 
tract both [O III] and H/3, assuming that the emission is 
70% of the [O III] flux (measured to within a factor of two, 
iTrager et al.ll2000a : lGrayes et al.ll200% . We then iterate these 
fits until the emission line flux has conyerged. In addition, we 
search for residuals around strong sky lines at 5200 and 5460 
A. 


From our iteratiye fits, and focused on the galaxy centers for 
simplicity, roughly two-thirds of the galaxies haye low-leyel 
emission detected, with a median EW of 0.2A, and a maximum 
of lA (calculated for those galaxies with detected H/3). The 
Gandalf measurements do not correlate yery strongly with our 
iteratiye fits. In the GANDALF fits, only half of the galaxies 
haye detections, with a median EW is 0.2A and a maximum H/3 
EW of 1.7A. As a means of quantifying our systematic errors, 
we rerun the stacking analysis described below on the Gandalf- 
subtracted spectra. As expected, only the stellar age changes 
significantly, being ~ 3 Gyr lower in the Gandalf stacks. But all 
other stellar population properties are yirtually identical (with 
the [Fe/H] shifting higher by a small amount to compensate the 
shift in stellar age). Thus, we present results based on the iter¬ 


atiye fits, but we caution that there is a rather large systematic 
uncertainty in the absolute stellar age. All other stellar popu¬ 
lation properties are robust to this modeling difficulty, and in 
general the relatiye ages are robust as well. We are currently 
working on more robust gas detection schemes using all lines 
in the spectra (V. Pandya et al. in preparation). 

We then use lick_ew (iGrayes & Schiayor] l2008h on the 
emission-line corrected spectr a. The indices are on a modi¬ 
fied Lick system presented by lSchiayonI (l2007l) based on flux- 
calibrated spectra. In order to demonstrate that we are on 
the same system, we compare the Lick indices from the flux- 
calibrated SDSS spectra (the inner 3") with those from the 
central 4f'2 fiber in our data. There is no net offset be¬ 
tween the two sets of indices in any case, with ((H/3EWs- 
H^EWms)/H/3EWms) =0.08±0.17, where S is SDSS and MS 
is the Mitchell Spectrograph. (Fe) and Mgb each haye a scatter 
of only ^ 10% and eyen smaller net offsets. 


4.3. Composite spectra 

While measuring Lick indices is a yery powerful technique 
for high S/N spectra, at the large radii that we are working, 
systematic effects such as small errors in sky subtraction and 
flux calibration can begin to cause large uncertainties in the 
Lick indices measured from indiyidual objects. Stacked spec¬ 
tra ayerage oyer sky subtraction and flux calibration errors in 
indiyidual systems, which occur at different wayelength s in 
each galaxy rest-frame (e.g., IGrayes et alJl2009l:lYanll2Ql Ih . Of 
course, yariations in stellar populations at a giyen tr* or mass 
are expected based on differences in accretion history (e.g., 
iHirschmann et alj|2015h and we are quite interested in these 
differences, particularly as a function of the dynamical prop¬ 
erties of the galaxies. We plan to implement full spectral fitting 
in the nea r future, which is more robust at low signal-to-noise 
ratio (e.g.. lChoi et al.ll2014i) . 

We know that st ellar population properties are a strong 
function of (e.g ., Worthey et alJ 1992t [Bender et al.lll993l : 
iTrager et al.ll2000a IGrayes et al.ll2009ll . Thus. we first diyide 
the galaxies into four stellar yelocity dispersion bins using our 
pPXF measurements to the central fiber. Of the 95 galaxies in 
our sample, we exclude eight. Most of these are from the low- 
mass sample and haye a bright star in the foreground making it 
difficult to reach a reasonable S/N. NGCl 167 in the MASSIVE 
sample is excluded because it may be a face-on disk galaxy. 
We will treat it more carefully in future work. The bins haye 
CT,< 220 km s“' (15 objects), 220 <(T*< 250 km s“* (21 ob¬ 
jects), 250 <(T,< 290 km s“' (33 objects), and all those higher 
than (T* > 290 km s“' (18 objects). Stacked spectra within 0.5Rc 
are shown in Figure|3] The resulting stellar population parame¬ 
ters are shown in Figure|4] We reach different physical radii for 
different bins as we run out of signal for the smaller and lower- 
mass galaxies. Also, the total number of galaxies included in 
the largest radial bin is typically ^ 30% smaller than the cen¬ 
ters, as yarious systematic effects such as foreground star con¬ 
tamination grow more seyere at low flux leyels. 

Below, we will also bin on stellar mass and group richness. 
While the detailed bins are different, the stacking technique de¬ 
scribed here is the same for these different sets of bins. 

To create the stack, we coadd the emission-line-subtracted 
spectra. We interpolate the rest-frame spectra onto a common 
wayelength grid. To ensure we know the dispersion of the fi¬ 
nal stacked spectrum, we then smooth each galaxy to a yalue 
that is 30% higher than the upper cr* limit of each bin. This 
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smoothing ensures that all galaxies go into the stack with the 
same effective dispersion. With this approach, we minimize 
small errors due to dispersion corrections to the Lick indices. 
An alternate procedure would be to smooth all templates, in 
all bins, to a high dispersion (e.g., 400 km s“'). We find very 
small (< 0.01 A) differences in the resulting indices if we adopt 
the latter approach, again confirming that our dispersion cor¬ 
rections are working (§4.2). When we make bins in Mk below, 
because of the wide range of tr, in each bin, we smooth all bins 
to 400 km s“'. 

We remove the continuum by dividing each spectrum by a 
heavily smoothed version of itself. This step simultaneously 
normalizes all spectra to the same level and ensures that differ¬ 
ences in continuum shape (whether real or due to small errors 
in sky subtraction or flux calibration) do not impact the final 
line strengths. We then calculate the median flux at each pixel, 
with rejection, alt hough we get very similar results using the bi¬ 
weight estimator (IBeers et al.|[l990h . We experiment with mul¬ 
tiplying the coadded spectrum by the median continuum before 
measuring indices, but the changes to the Lick indices are neg¬ 
ligible. 


4.4. Uncertainties 

To determine the level of variation in the composite spectra, 
we generate 100 boot-strapped composite spectra by randomly 
drawing from the total list of galaxies in that bin, with replace¬ 
ment. We measure Lick indices from each of these 100 trial 
spectra. We then assign errors on the Lick indices that enclose 
68 % of the Lick indices measured from the 100 trials. There¬ 
fore, the size of the error is most directly related to the variance 
in parameters over the population in that bin, rather than mea¬ 
surement or modeling uncertainty. The systematic uncertainties 
are not shown. Age measurements have large systematic uncer¬ 
tainties due to emission infill (§4.2), while nitrogen and calcium 
are particularly uncertain due to their dependence on blue spec¬ 
tral features and assumed carbon and oxygen abundances. 

In general, we report the measurement from the primary 
stacked spectrum, and the errors derived from the boot-strapped 
spectra. There are a few cases where the Lick index measured 
from the primary stacked spectrum does not fall within the stel¬ 
lar population grids (usually because the H/3 index is slightly 
too low). In these cases, we use the median index value from 
the 100 trials as the final answer. These cases are indicated with 
open symbols. 

5. RADIAL VARIATIONS IN STELLAR POPULATIONS 
5.1. Bins of a ^ 

In Figure |4] we show radial trends in the measured age, 
[Fe/H], and abundance ratios as a function of physical (left) 
and /?e-scaled (right) radii. Our default models make the as¬ 
sumption that [O/Fe] is enhanced like [Mg/Fe], since they are 
both a elements. However, in Figure|4]we indicate with dashed 
lines the C, N, and Ca abundances that result for an assumed 
solar [O/Fe]. While the zeropoints of [C/Fe] and [N/Fe] both 
decline, the radial trends will not change unless [O/Fe] (unlike 
[Mg/Fe]) changes with radius. 

We fit a power-law relation between radius and each stel¬ 
lar population property. We anchor the relation at the cen¬ 
ter of our radial coverage. Eventually, when we have reli¬ 
able //-band light profiles for each galaxy, we will calculate 
mass-weighted stellar population properties, but at present, for 


a given stellar population property we fit a log-linear rela¬ 
tionship : X = A log(/?//? 3 )-l-B, where is the third bin, cor¬ 
responding either to 3-6 kpc or 1 - l.SRg. The quantity A rep¬ 
resents the radial gradient per log radius. The quantity B repre¬ 
sents an effective stellar population property at the center of our 
radial coverage. The fits to A and B for each relation are shown 
in Tables 1 & 2. 

The strongest radial gradients are found between [Fe/H] and 
[C/Fe], which both decline with radius. Age shows a decline 
in the highest dispersion bin, but because that effect is weaker 
when we use the alternate emission line correction, we treat age 
gradients with extra caution. All other abundance ratios mea¬ 
sured here are consistent with remaining flat over the full radial 
range, aside from a 2 (t increase in [Ca/Fe] at large radius that is 
only seen in the largest cr* bin. 

We mrn to trends between stellar populations and cr*. Our 
central bins are shown as open circles in Figure |3 We re¬ 
cover well-known trends between stellar population properties 
and stellar velocity dispersion for galaxy centers. Galaxies 
with higher stellar velocity dispersions ha ve older stellar ages 


Worthev 

2004; Thomas et alJ 

20051; 

Sanchez-Blazauez et al.l 

20061: Graves et al. 2007: Smith et alJ 

20091; Price et all 2011; 

Johansson et al.l2012; Worthev 

et al.l2014; Conrov et alj2014l). 


One interesting exception is [N/Fe], where we see a decline 
with Q-j,. This tre nd is at odds with our previous finding 
(iGreene et al.ll2013[), as well as most other work on the topic 
(although see also Kelson et al.ll2006l) . More detailed work is 


needed to confirm this trend, particularly given possible uncer¬ 
tainties in flux calibration in the blue for the Mitchell spectra. 
Recall also that the very old ages found in the very centers of 
these galaxies have 2-3 Gyr error bars, and thus are consistent 
with the age of the universe. 

Thanks to our spatial coverage, we can go beyond the trends 
between cr* and stellar populations in the galaxy center, and 
look at how these trends evolve as we look to larger radius, 
both in physical and R(,-scaled units. The crosses in Figure |5] 
show the trends between cr* and stellar population properties 
measured at 3-6 kpc. In general, the same trends are seen with 
cr* in the central and 3-6 kpc bins ([N/Fe] is again an excep¬ 
tion). Fit values are included in Table 3. [C/Fe] shows the most 
significant correlation with cr*, with a positive slope detected 
at 8cr significance. We find a weak trend with [cr/Fe] (3.5cr) 
and a strong trend with age (7 cr; not shown). The measured 
slopes between cr* and [Fe/H], [Ca/Fe], and [N/Fe] are all con¬ 
sistent with zero; that is, the effective value of these parameters 
is similar for all cr* bins. 

We can also evaluate trends between cr* and stellar pop¬ 
ulation properties evaluated beyond R^ (taken here to be the 
1 - L5Re bin). In this case, the physical radii increase for the 
higher-dispersion bins. Interestingly, we find that the trends 
with cr* are weaker when evaluated as a function of Re- Only 
stellar age (7 cr) and [C/Fe] (much more weakly at 4cr) correlate 
positively with cr*. Apparently, trends between cr* and abun¬ 
dance ratios are strongest when measured within small phys¬ 
ical radii. If galaxies indeed form a compact core rapidly at 
high redshift, then we e xpect this inner region (< 5 kpc; e.g., 
Ivan der Wei et al.ll2014l) to depend most strongly on cr*. 

Before we interpret these observed trends in more detail, we 
examine differences between bins in cr* and stellar mass. 
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Fig. 5.— Left'. We compare the relation measured from the coadded central fibers (open circles; dashed lines) with those from fits to the profiles weighted to 
1 — l.5Re (filled circles; solid lines) and 3 — 6 kpc (crosses; long-dashed lines). We plot the fitted zeropoints (B) for each chemical property fitted versus radius 
(Tables 1 & 2) as a function of cj*. [Z/H] is estimated as [Z/H]=[Fe/H]+0.94[Q;/Fe], which assumes that [0/Fe] tracks [Mg/Fe]. The power-law fits shown here have 
the form X = A log (cr*/250km s“^) + B for each stellar population property X (Table 3). We see that the central and 3-6 kpc fits are generally similar and show 
the strongest trends with cr*, while when we scale to Re the trends are washed out, supporting a scenario where the inner 5 kpc collapsed quickly and scales with 
the galaxy potential, while the outer parts are accreted from lower-mass systems. Right. As above, but fitting X = A log {Mk/ — 253 mag) + B. Note that the central 
measurement is missing for the highest Mk bin because the H/3 measurement fell off the grid. In general the dispersion in galaxy properties is lai'ger in fixed Mk 
bins, and the trends weaker. The largest exception is [C/Fe] which shows a stronger correlation with Mk than cr*. 
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Table 1. Gradient Fits (Re) 


Stack 

Aage 

Bage 

A.[Fe/H] 

B[Fe/H] 

A[„/Fe| 

B[a/Pe] 

A[c/Fb] 

B[C/Fe] 

■A-[N/Fe] 

B[N/Fe] 

A-[Ca/Fe) 

B[Ca/Fe] 

<^<220 

-0.36±0.33 

7.41±0.82 

-0.39±0.31 

-0.37±0.08 

0.15±0.12 

0.32±0.03 

-0.23±0.08 

0.12±0.02 

-0.21±0.12 

0.80±0.04 

0.10±0.13 

0.12±0.04 

<^220-250 

-0.03±0.17 

10.2±0.48 

-0.34±0.10 

-0.32±0.03 

0.01±0.06 

0.29±0.02 

-0.34±0.07 

0.18±0.02 

-0.16±0.09 

0.71±0.03 

0.06±0.06 

0.08±0.02 

<7250-290 

-0.34±0.21 

10.8±0.77 

-0.28±0.18 

-0.34±0.05 

-0.04±0.10 

0.31±0.03 

-0.43±0.07 

0.11±0.02 

0.13±0.14 

0.83±0.04 

0.03±0.08 

0.07±0.03 

<7290-400 

-0.35±0.11 

12.9±0.34 

-0.46±0.11 

-0.40±0.02 

0.05±0.06 

0.35±0.02 

-0.32±0.06 

0.22±0.01 

0.09±0.10 

0.78±0.03 

0.17±0.08 

0.14±0.03 

Afx>-24.5 

-0.37±0.24 

9.59±0.56 

-0.27±0.11 

-0.35±0.03 

0.04±0.05 

0.34±0.02 

-0.33±0.06 

0.13±0.01 

-0.21±0.10 

0.86±0.03 

0.02±0.05 

0.07±0.02 

Afx-24.5:-25.5 

-0.09±0.26 

11.3±0.74 

-0.44±0.23 

-0.32±0.05 

-0.04±0.11 

0.26±0.03 

-0.36±0.06 

0.16±0.01 

-0.27±0.11 

0.76±0.03 

-0.01±0.11 

0.04±0.02 

Afx-25.5:-25.8 

-0.34±0.17 

9.01±0.56 

-0.37±0.13 

-0.30±0.04 

0.03±0.08 

0.29±0.03 

-0.31±0.08 

0.17±0.02 

-0.03±0.12 

0.72±0.04 

0.09±0.10 

0.09±0.03 

Afx>-25.8 

-0.29±0.29 

13.0±0.89 

-0.54±0.17 

-0.52±0.06 

0.19±0.08 

0.43±0.03 

-0.36±0.09 

0.20±0.03 

0.43±0.12 

0.85±0.04 

0.10±0.10 

0.17±0.03 

Mh 

-0.17±0.21 

8.85±0.68 

-0.51±0.16 

-0.43±0.05 

0.19±0.08 

0.41±0.03 

-0.24±0.07 

0.20±0.02 

0.16±0.09 

0.84±0.02 

0.27±0.06 

0.14±0.02 

M- 

-0.24±0.14 

8.93±0.38 

-0.30±0.08 

-0.28±0.02 

-0.06±0.04 

0.28±0.01 

-0.37±0.04 

0.17±0.01 

-0.23±0.07 

0.85±0.02 

-0.003±0.03 

0.003±0.01 

ENVfc 

-0.23±0.27 

9.78±0.83 

-0.44±0.14 

-0.42±0.04 

0.10±0.08 

0.38±0.03 

-0.31±0.06 

0.19±0.02 

0.06±0.10 

0.87±0.03 

0.19±0.07 

0.10±0.02 

ENVj 

-0.29±0.14 

7.67±0.43 

-0.21±0.10 

-0.26±0.03 

-0.09±0.06 

0.27±0.02 

-0.38±0.05 

0.14±0.01 

-0.13±0.11 

0.82±0.03 

-0.04±0.07 

-0.002±0.02 

BCGh 

-0.31±0.07 

9.40±0.03 

-0.40±0.20 

-0.41±0.06 

0.15±0.07 

0.37±0.02 

-0.31±0.07 

0.18±0.03 

0.08±0.12 

0.84±0.03 

0.19±0.06 

0.10±0.02 

BCGi 

-0.16±0.18 

8.49±0.57 

-0.40±0.11 

-0.34±0.04 

0.06±0.07 

0.35±0.02 

-0.31±0.07 

0.16±0.03 

0.04±0.10 

0.90±0.02 

0.11±0.06 

0.06±0.02 


Note. — Fits to stellar population paramaters as a function of radius in units of Re- For each parameter X, we fit a relation X = A log(i?/1.5ile) + B. l.Sile is the third radial bin, so B 
represents an effective value for the stellar population property, while A is the radial gradient. These fits are done to stellar population properties derived from bins of (in order from top to 
bottom) <7*, Mji, Mji at fixed <7*, environment or halo mass at fixed <7*, and then halo mass for only the brightest galaxy in the halo. 


Table 2. Gradient Fits (Physical) 


Stack 

Aage 

Bage 

A[Fe/H| 

B|Fe/H| 

A|„/Fe| 

B|a/Fe] 

A[C/Fe] 

B[C/Fte] 

A[N/Fel 

B[N/Fel 

A[C./Fe] 

B[Ca/Ffe] 

<7<220 

-0.28±0.24 

7.75±0.83 

-0.20±0.10 

-0.31±0.04 

-0.02±0.03 

0.27±0.02 

-0.27±0.04 

0.11±0.02 

-0.23±0.07 

0.81±0.03 

0.03±0.05 

0.08±0.03 

<7220-250 

-0.06±0.12 

9.83±0.48 

-0.27±0.06 

-0.29±0.02 

0.00±0.05 

0.28±0.02 

-0.19±0.04 

0.21±0.02 

-0.24±0.08 

0.69±0.04 

-0.04±0.05 

0.04±0.02 

<7250-290 

-0.42±0.11 

11.7±0.34 

-0.10±0.04 

-0.29±0.02 

-0.02±0.03 

0.30±0.01 

-0.24±0.03 

0.21±0.01 

0.04±0.07 

0.79±0.02 

0.05±0.03 

0.08±0.01 

<7290-400 

-0.47±0.12 

13.0±0.43 

-0.19±0.05 

-0.28±0.02 

-0.02±0.03 

0.34±0.02 

-0.25±0.03 

0.29±0.01 

-0.02±0.06 

0.73±0.02 

0.05±0.03 

0.08±0.02 

Mx>-24.5 

-0.42±0.16 

9.78±0.64 

-0.20±0.07 

-0.36±0.03 

-0.02±0.03 

0.33±0.02 

-0.29±0.04 

O.lliO.Ol 

-0.09±0.05 

0.88±0.03 

0.04±0.03 

0.08±0.02 

A^X-24.5:-25.5 

-0.11±0.12 

11.2±0.58 

-0.26±0.15 

-0.32±0.05 

-0.01±0.08 

0.28±0.03 

-0.25±0.04 

0.20±0.02 

-0.19±0.10 

0.77±0.04 

0.02±0.08 

0.05±0.03 

A^X-25.5;-25.8 

-0.14±0.13 

9.74±0.42 

-0.17±0.06 

-0.22±0.02 

-0.04±0.03 

0.28±0.01 

-0.25±0.04 

0.22±0.01 

-0.09±0.06 

0.73±0.02 

0.03±0.04 

0.10±0.02 

A^X>-25.8 

-0.76±0.23 

13.1±0.49 

-0.22±0.12 

-0.28±0.03 

0.08±0.08 

0.35±0.02 

-0.27±0.05 

0.30±0.01 

0.35±0.14 

0.70±0.03 

0.21±0.10 

0.12±0.02 


-0.17±0.11 

9.06±0.35 

-0.15±0.09 

-0.26±0.04 

0.00±0.04 

0.34±0.02 

-0.22±0.03 

0.26±0.01 

-0.10±0.06 

0.78±0.02 

0.04±0.05 

0.06±0.02 

m; 

-0.24±0.10 

9.37±0.39 

-0.24±0.04 

-0.30±0.02 

-0.02±0.02 

0.28±0.01 

-0.25±0.03 

0.18±0.01 

-0.14±0.05 

0.87±0.03 

0.04±0.02 

0.03±0.01 

ENVft 

-0.27±0.13 

10.1±0.46 

-0.17±0.09 

-0.32±0.04 

-0.01±0.04 

0.33±0.02 

-0.21±0.03 

0.25±0.01 

-0.05±0.05 

0.83±0.02 

0.07±0.05 

0.06±0.02 

ENVj 

-0.07±0.11 

8.51±0.38 

-0.20±0.06 

-0.25±0.02 

-0.03±0.03 

0.28±0.02 

-0.26±0.07 

0.15±0.02 

-0.15±0.10 

0.86±0.04 

-0.04±0.08 

-0.01±0.03 

BCG-h 

-0.25±0.04 

10.5±0.02 

-0.18±0.08 

-0.34±0.04 

-0.01±0.04 

0.33±0.02 

-0.20±0.03 

0.25±0.01 

-0.05±0.07 

0.80±0.02 

0.04±0.05 

0.04±0.02 

BCG-1 

-0.11±0.09 

8.38±0.31 

-0.20±0.06 

-0.27±0.02 

-0.00±0.04 

0.33±0.02 

-0.21±0.02 

0.22±0.02 

-0.12±0.06 

0.86±0.03 

0.01±0.05 

0.02±0.02 


Note. — Fits to stellar population paramaters as a function of physical radius. For each parameter X, we fit a relation X = A log (H/i? 3 _ 6 kpc) "b B. fls-ekpc is the third radial bin, so B 
represents an effective value for the stellar population property, while A is the radial gradient. These fits are done to stellar population properties derived from bins of (in order from top to 
bottom) (7*, Mx) Mx at fixed cr*, environment or halo mass at fixed cr,, and then halo mass for only the brightest galaxy in the halo. 
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Fig. 6.— Radial gradients in age, [Fe/Fl], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] as above, but now in bins of Mj( (mag; a proxy for stellar mass) rather than stellar 
velocity dispersion. Open symbols indicate cases where the primary measurement fell off the stellar population grid and we have used the median of the distribution 
of boot-strap stacks instead. The inversion program failed to converge for the [N/Fe] and [Ca/Fe] measurements for the outermost high-Mx Rc-scaled point, so 
no measurement is shown. As indicated in the key, the bins have Mk > —24.5 (squares), —24.5 > Mk > -25.5 (pentagons), —25.5 > Mk > -25.8 (triangles), and 
Mk < —25.8 mag (circles) respectively. We fit the radial gradients with a power law of the formX = A log/R/Rsf-rB for each stellar population parameter X, where 
R} is either 3 — 6 kpc or 1 — 1.5Re. Fits to the most massive (long-dashed lines) and least massive (dotted lines) galaxies are shown here; all fitted coefficients are in 
Tables 1 & 2. 
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5.2. Bins ofM* 


There is now a preponderance of evidence that stellar pop- 
ulati on properties correlate most strongly with at their cen¬ 

ters (iBender et al.ll993tlTYager et alj2000btlCiraves et al.l2009t 


IWake et ^l2012h . However, there is little work comparing ra - 

dial trends as a function of mass and tr* (ISnolaor et al.ir2010r) . 
We therefore create f our bins of A'-band magnitude, based o n 
2MASS photometry (iJarrett et al.ll2003l : ISkrutskie et alJl2006l) : 


Mic > -24.5 (21 galaxies), -24.5 > Mk > -25.5 (23), -25.5 > 
Mfc > -25.8 (24), andM/f < -25.8 mag (19; Figure|6ll. We note 
that there is a wide range of stellar velocity dispersion at each 
mass bin (Figure |2|, so this binning scheme is truly different 
from the dispersion bins presented above. 

As above, we fit each stellar population property as a func¬ 
tion of radius, tabulated in Tables 1 & 2. Each effective stellar 
population parameter is shown as a function of Mk in Figure |3 
(right). Starting as in §5.1 by focusing on the central measure¬ 
ments, we see that in general, there is more dispersion in stellar 
population properties (particularly metallicity and [Mg/Fe]) in a 
given M/f bin as compared with a given cr* bin. In particular, the 
second Mk bin (-24.5 to -25.5 mag) shows considerably larger 
variance than the other bins. We believe that this larger error bar 
reflects a genuine increase in the spread in stellar populations 
in this bin. At yet lower Mk, we no longer have a representative 
sample of objects. We also generally see weaker trends in the 
stellar populations with Mk than with tr*, particularly in phys¬ 
ically scaled bins (Table 3). The exception is [C/Fe], which 
shows amore significant correlation with M/f (12cr). We see no 
correlation with age or [a/Fe]. There is, however, a marginal 
correlation with [Fe/H] (3 a) and a corresponding weak corre¬ 
lation with [Z/H] (also 3 a). If we instead consider measure¬ 
ments from l.SRe, we And no significant trends between stellar 
population properties and Mk (not even stellar age). 


5.2.1. Bins in M* at fixed cr* 

We now ask whether, at fixed cr*, there are residual trends 
in the stellar population gradients as a function of stellar mass 
(Figure^. We create two luminosity bins divided at Mk = -25 
mag (which approximately divides the sample in two). We then 
enforce a matched distribution in cr* between 200-330 km s“^ 
by creating 100 stacks for each mass bin, each having an identi¬ 
cal distribution of cr* and each including a total of 30 galaxies. 

We And small but detectable differences in the two mass bins. 
Specifically, while the average stellar population ages are both 

9.5 Gyr, we And that the high-mass bin is more a-enhanced 
by 0.1 dex with (possibly) correspondingly lower [Fe/H] at 
Axed Re (Table 1). The C and N abundances and radial proflles 
are consistent with each other, but there is a hint that the more 
massive galaxies are also more Ca-enhanced at Axed cr*. All of 
these trends could be qualitatively explained if the timescale for 
star formation were shorter in the more massive systems, such 
that [Fe/H] is lower while [a/Fe] is higher and [Ca/H] follows 
[Fe/H]. 

The other possible explanation is that at Axed cr*, more mas¬ 
sive galaxies are physically larger. This would explain why at 
a Axed fraction of R e , the more massive galaxies have lower 
[Fe/H]. [Graves et al.l (l2009h And that galaxies with the lowest 
central surface brightness also have the oldest ages and lowest 
[Fe/H] values, a similar trend to that seen here, albeit only for 
galaxy centers. When we have robust size measurements for 
our sample galaxies we will revisit the question of stellar pop¬ 
ulations through the Fundamental Plane, also including radial 


stellar population information. 

5.3. Interpreting Stellar Populations Using Radial Trends 

There are a number of interesting trends seen in the stellar 
populations of elliptical galaxy centers that challenge our un¬ 
derstanding, particularly when compared with Galactic trends. 
It is our hope that adding radial information will shed new light 
not only on the assembly history of ellipticals, but also on the 
nucleosynthetic yields that lead to these observed trends. 

There are a few important caveats to keep in mind as we in¬ 
terpret the observations. First, our stellar ages, particularly at 
the center, are subject to large (^ 2-3 Gyr) uncertainties. Thus, 
we caution against over-interpreting the radial trends in age at 
present. Some work covering l arger dynamic range in r adius do 
detect clear age gradients (e.g.. lLa Barbera et al.ll2012l) . but we 
do not go to large enough radius to detect the very low metal¬ 
licity true h alo component tha t has been seen in a f ew nearby 
cases (e.g., [Harris et aDl2007t I Williams et al.l l20 ill) . Second, 
we remain cautious about the nitrogen measurements given the 
possible flux calibration difflculties at the blue end of the spec¬ 
trum. Third, recall that we are examining average trends in 
stacked spectra. Undoubtably there are interesting exceptions 
to all of these trends, and we plan to study the full range of 
parameters in future work. 

We revisit our expectations for stellar population trends as 
a function of radius in light of the recent picture that galaxies 
form in two “phases”, an initial burst of in-situ star formation 
creating the central component, followed by late- time accre¬ 
tion o f smaller, fluffler units at larger radius (e.g., iNaab et ^ 
l2009ll . For galaxies in our mass range, the typical sizes of the 
central components are measured to be 2-5 kpc at z 2 (e.g., 
Ivan der Wei et al.l 1201 4ll . When we look at the stellar popu¬ 
lations at this typical radius (Figure |5]l, we And that age and 
[a/Fe] depend on cr* rather than Mk- Since we expect galax¬ 
ies that form early to be denser and have higher cr*, a stronger 
correlation between age and star formation timescale with cr* 
seems natural. In contrast, [Fe/H], [C/Fe], and [N/Fe] depend 
more on the total stellar mass of the system. 

In addition to considering a Axed physical radius, we look 
at stellar population trends at ^ l.SRe (Figure |5]l. According 
to simulations, more massive galaxies are increasi ngly domi¬ 
nated by accreted stars at radii beyond ^ 5 kpc (e.g.. lOser et al.l 
I 2 OIOI) . Therefore we might expect the correlation between cr* 
and stellar population properties to decrease when taken over 
the bulk of the stellar population. This is what we observe. 
Beyond ^ Re we And no strong trend between cr* or Mk and 
abund ances or abundance ratios. This result builds on what we 
saw in iGreene et akl (120131) . There, we emphasized that the stel¬ 
lar populations beyond Re in massive galaxies tend to have low 
[Fe/H]-0.5 dex and high [a/Fe]^ 0.3 dex, stellar populations 
that are not seen in the centers of any galaxies today (see also 
iBender et al.ll2015T) . Here we see that beyond Re, galaxies over 
a relatively wide range in cr* and Mk have similar abundances 
and abundance ratios, as expected if the more massive galaxies 
were built by accreting the less massive. Only stellar age is still 
seen to increase at higher ct* when examined at ^ Re', we await 
better age measurements to verify this result. 


Another ongoing discussion in the literature regards carbon. 
As is seen here, [C/Fe] is observed to increase with cr* in ellip- 


tical galaxy centers (e. 

2 .. Tracer et alJl998HGraves et alJl2007l: 

iJohansson et al.ll2012l 

Conroy et al.ll2014l). This is in contrast 


to the behavior of carbon in the Milky Way or Local Group 
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Table 3. Galaxy Correlations at Varying Radii 


Radius 

A[Fe/H] 

B[Fe/H] 

A[z/h] 

B[z/h] 

A[„/Fe] 

B[„/Fel 

A[C/Fe] 

B[C/Pe] 

A[N/Fe] 

B[N/Fe] 

A[Ca/Fe] 

B[Ca/Fe] 

(T* Fit 

Center 

1 - 1.5Re 

3 — 6 kpc 

-0.23±0.43 

-0.56±0.28 

0.13±0.20 

-0.17±0.03 

-0.34±0.02 

-0.29±0.01 

0.28±0.47 

0.02±0.36 

0.63±0.22 

0.12±0.03 

-0.06±0.03 

-0.02±0.02 

0.44±0.15 

0.34±0.18 

0.41±0.12 

0.31±0.01 

0.32±0.01 

0.30±0.01 

0.93±0.22 

0.54±0.11 

0.96±0.11 

0.37±0.01 

0.16±0.01 

0.19±0.01 

-1.25±0.46 

0.13±0.26 

-0.33±0.20 

0.83±0.03 

0.77±0.02 

0.77±0.01 

-0.05±0.22 

0.19±0.25 

0.15±0.18 

0.06±0.01 

0.09±0.01 

0.07±0.01 

Mk Fit 
Center 

1 - 1.5Re 

3 — 6 kpc 

-0.04±0.04 

0.02±0.03 

-0.07±0.02 

-0.15±0.04 

-0.36±0.02 

-0.27±0.01 

-0.03±0.05 

0.01±0.03 

-0.06±0.02 

0.15±0.05 

-0.05±0.03 

0 .01±0.02 

0 .01±0.02 

-0.02±0.01 

O.OliO.Ol 

0.32±0.02 

0.34±0.01 

0.31±0.01 

-0.07±0.03 

-0.03±0.01 

-0.09±0.01 

0.38±0.02 

0.17±0.01 

0 .22±0.01 

0.10±0.04 

0.03±0.02 

0.09±0.02 

0.S0±0.04 

0.79±0.02 

0.76±0.01 

-0.02±0.01 

-0.04±0.02 

-0.02±0.01 

0.07±0.02 

0.09±0.01 

O.lOiO.Ol 


Note. — Fits to stellar population paramaters as a function of it* or Mjr, weighted to different radii. For each stellar population parameter X, we fit a relation X = A log ((T*/250 kms~^} + B 
(top) and X = A log (MkI — 25.3mag) B (bottom). X comes from the central binned data, from the fit in Table 1 weighted to 1 — 1.5i?e, and finally from the fits in Table 2, weighted to 
3 — 6 kpc. 


dwarf galaxies (e.g. jKirbv et al.ll2015h . To get such super-solar 
carbon levels in the short timescales implied by the high ratio 
of [a/Ee], carbon must come not only fro m intermediate-mass 
(AGB) stars but also from massive stars (iGraves et all 1200'^ 
iTang et aTll2014l) . Carbon yields from massive stars are thought 
to increase with incre asing metallici ty due to increased mass- 
loss from winds (e.g.. lMaedeijll992h . These same high yields 
at high metallicity are als o invoked to explai n abundance trends 
in Milky Way stars (e.g. jHenrv et al.lbOOOh . If true, we would 
expect that as the metallicity decreases outwards in these el¬ 
liptical gala xies, the IC/Ee] wo uld also decrease. This is what 
we observe dGreene et alj|2013l) . Interestingly, we find that the 
[C/Ee] gradient follows the decline in [Ee/H] in both the cr* and 
Mk bins. 

There is also considerable debate in the literature about the 
origin and behavio r of nitrogen in elliptica l galaxy centers (e.g., 
iKelson et alJuPOa IJohansson et alJl^l^ . Our surprising re¬ 
sults are two-fold. Eirst, [N/Ee] is remarkably super-solar. Even 
assuming solar [O/Ee] (at odds with the observed [Mg/Ee]) 
[N/Ee] is found to be three times the solar value. Second, since 
the N is produced by C through the CNO cycle, their differ¬ 
ent behavior with radius is non-intuitive (and may point to a 
changing O abundance as well). As discussed above, the flux 
calibration at the blue end of the spectrum, containing CN, is 
quite uncertain. We do And good agreement between CN as 
measured from the SDSS spectra and the central Mitchell fiber. 
We also confirm that the CN measurements do not depend on 
how we treat the continuum. Including or excluding the overall 
continuum level in the stacks changes the CN 1 measurement by 
< 0.005 mag, resulting in very small < 0.02 dex in the [N/Ee] 
ratio. To really confirm these high N abundance ratios at large 
radius, we would like to perform full spectral modeling to miti¬ 
gate the impacts of blending and the uncertainties introduced by 
oxygen. In the meantime, it is interesting to note other stellar 
systems that display very super-solar nitrogen abundance ra¬ 
tios. Eor instance, in globular clusters, the wide range in [N/Ee] 
strongly suggests p re-enrichment by a previous early epoch of 
star formation (e.g.. lCohen et al.l[2005l) . 


Einally, we come to Ca. While nominally an a-element, it 
has long been known, based on both Ca4227 in the blue and 
the calcium triplet (CaT) index in the red, that Ca is under- 


abundant with respect to the other a elements (e.g.. Cohen 

1979; 

IVazdekis et al. 

1997;|Worthev 1998 

;|Proctor & Sansom 

2002; 

Terlevich et al. 

1990t Peletier et al. 

Il999l Sagliaetal. 

2002; 

Thomas et al.l 

20031 IChoi et al. 2014. Like FEe/Hl. 


[Ca/Ee] shows no dependence on cr*. There are a number of ex¬ 
planations in the literature for the CaT mea surements, including 
changes in the initial mass function (e.g.. ICenarro et al]l2004ll 


or a minority metal-poor population (ISaglia et al.ll2002h . How¬ 
ever, to explain both the blue and red index behavior, it is more 
natural to presume that Ca behaves like an Ee-peak element be- 
cause it is predominantly produced in Type la supernoy ae (e.g., 
IWorthevll998l : IWorthev et al.l201 itl^nrov et alJ2014l) . As ex¬ 
pected in that case, we measure a flat [Ca/Ee] ratio with radius. 
The one intriguing difference is the possible increase in [Ca/Ee] 
at large radius, which is worth pursuing. 

5.4. Bins of Group Richness 

While there are well-documented differences in the mor- 
phol ogical mix of g alaxies as a function of local galaxy den¬ 
sity dDressleJ 119801) . the observations of environmental dif- 
ferences in stellar population proper t ies are quite subtle (e.g ., 
|Th omas et a n i200-i IZhu eTZl bcilOt lEackner & Gunnll20li . 
Early stu dies found evidence for yo unger ages in ‘field’ galax¬ 
ies (e.g.. iTerlevich & Eorbesll2002l) . A number of other stud¬ 
ies report no change in scaling relations between stellar pop- 
ulati on parameters and (7„ as a function of local environment 
(e.g.. lKuntschner et alJl2002l : [Bernardi et alJl2006l) . but do And 
a larger fraction of “rejuvenated” galaxie s with recent star for¬ 
mation in low-densi ty environments (e.g.. lAnnibali et aT]l2007l : 
iThomas et alj|201^ . Recently, the samples have grown large 
enough to evaluate not just the average properties of field and 
cluster galaxies, bu t to control for stellar and halo mass (e.g., 
iPasQuali et alJ201^ . Eor instance, Pasquali et al. And that satel¬ 
lite galaxies at fixed mass grow older and more metal rich as 
their host halo mass increases. 

There is precious little literature on radial gradients in stel¬ 
lar populations as a function of group richness, although a few 
photometric studies A nd steeper metallicity gradients in lower - 
density environments (iKo & Imll2005l : iLa Barbera et alJl2005h . 
Give n our large samp le and wide range of measured group rich¬ 
ness (iMa et al.ll2014l) . we are in a unique position to examine 
radial stellar population trends with environment. 

To divide the samp le by richness, we use the group catalog 
of ICrook et alJ (l2007h . When we divide our galaxies based on 
group richness alone, the distributions in cr* do not match. In¬ 
stead, there are more high cr, galaxies in richer groups. Since 
CT, is strongly correlated with stellar population properties, we 
must match the cr* distributions across different halo mass bins. 
Therefore, we divide the sample in half based on the number of 
neighbors; “low” comprises galaxies with no more than three 
companions of > L* while “high” comprises the rest. The raw 
distribution of tr* have median (tr*) = 240 km s“' for the low- 
density bin and (cr*) = 260 km s“' for the high-density bin. 
Since there are not sufficient numbers to further subdivide the 
galaxies into bins of cr*, we create 100 stacks, drawing from 
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Fig. 7.— Radial gradients in age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] as above, but now in cr*-matched bins of Mf( (mag; a proxy for stellar mass). 
Filled symbols have Mk < -25 mag while open symbols have Mk > —25 mag. We fit the radial gradients with a power law of the form X = A log(R/R 3 ) + B for 
each stellai' population parameter X, where is either 3-6 kpc or 1 - l.5Re. The fits to the high-mass (long-dashed lines) and low-mass (dotted lines) bins are 
shown here and in Tables 1 & 2. 
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Fig. 8.— Radial gradients in age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] as above, but now in bins of group richness, controlling for the distribution of cr*. 
Low density (open squares) have three or fewer neighbors with L> L*, while filled circles include everything else. We fit the radial gradients with a power law of 
the form X = A log(/?//? 3 ) + B for each stellar population parameter X, where is either 3-6 kpc or 1 — l.5Re. The fits to the rich (long-dashed lines) and poor 
(dotted lines) bins are shown here and in Tables 1 & 2. 
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the objects with cr* between 200 and 330 km s“'. We force the 
low-density and high-density stacks to have the same number of 
objects (24 in this case) with the same distribution of cr*. We do 
not have sufficient numbers to match on as well; two-thirds 
of the rich galaxies are also in the brighter half of the sample. 
The results are shown in Figure[8] 

There are only very slight differences between the two galaxy 
samples divided by group richness. The high-density stack is 
slightly older, has slightly lower [Fe/H], and slightly higher 
[g/Fel, similar to some previous studies (e.g ., iBernardi et al.l 
l2006l:lQemens et al.l2009tlCoot)er et al.l201ol) . The basic inter¬ 
pretation is that objects found in the highest density peaks today 
likely formed earlier, and thus have older ages, higher [a/Fe], 
and slightly lower [Fe/H] (although they maintain roughly solar 
metallicity overall). Again, the other possibility is that at fixed 
tr*, galaxies in denser environments tend to be slightly larger. 
We plan to control for galaxy size in future work. Turning to the 
other light elements, we find that [C/Fe] is marginally higher in 
the high-density bin, while [N/Fe] and [Ca/Fe] are comparable 
between the two. Finally, it is interesting to note that the sample 
variance is larger for the low-density sample. 

We see marginal evidence for a steeper [Fe/H] gradient in the 
higher halo-mass bin as a function of effective radius (FigurejS) 
right; -0.44 ± 0.14 for the higher halo mass bin, -0.2 ± 0.1 for 
the lower halo mass bin). The slope difference is only signifi¬ 
cant at ^ 2 cr. [of/Fe] also shows a very marginal difference in 
the opposite direction, with slopes of 0.1 ±0.08 and -0.09 ± 
0.06 for the higher and lower halo mass bins respectively. Be¬ 
cause of the anti-correlation between [Fe/H] and [cr/Fe], the 
resulting gradients in [Z/H] for the two bins are comparable 
(-0.3 ±0.1 in both cases). If this result is confirmed, it perhaps 
suggests that at a given cr», galaxies in richer environments are 
more compact, and thus show steeper [Fe/H] gradients. While 
we do see marginal differences between [Fe/H] and [a/Fe], we 
do not confirm results from previous photometric surveys that 
reporte d steeper metal l icity gradients in lower density environ- 
m ents dK o & Imll 2()()5 : La Barbera et alJl200^ . 

iLa Barbera et al.l (l2014t) argue that in addition to looking for 
trends as a function of halo mass, we should also be dividing the 
samples into “central” (the most massive galaxy in a halo) and 
“satellite” galaxies. We therefore repeat the stacking exercise, 
but this time removing all satellite galaxies (Figure |9]l; we do 
not have sufficient numbers of satellite galaxies to stack them 
alone. Interestingly, when we examine central galaxies alone, 
the marginal stellar population differences discussed above van¬ 
ish (although we note that our statistical power is lessened by 
the smaller sample size as well). Perhaps the differences be¬ 
tween the two samples are driven by differing fractions of satel¬ 
lite galaxies in the two environmental stacks, but better statistics 
are needed before we can be sure. On the other hand, studies 
of individual brightest cluster galaxies in rich clusters find a 
wide spread in central proper ties such as age and [a/Fe] (e.g., 
lOliva-Altamirano et al.ll20f5l), as well as evidence for distinct 
accretion episodes (e.g., Coccato et al.ll2010l 1201 ih . In future 
work we will investigate in more detail the spread in central 
galaxy properties as a function of halo mass. 

We do not reproduce the trend found by La Barbera et al. that 
the central galaxies are younger in the larger halos. Instead, 
when we focus on central galaxies only, we see no significant 
difference between the stellar populations of the two groups at 
the 0.1 dex level. On the other hand, “low-mass” halos in the 
La Barbera study have M* < 10'^ ^ Mq. This is considerably 


lower than the likely halo mass of our galaxies, given that their 
stellar masses reach M* ss lO'^ Mq. Furthermore, the differ¬ 
ences reported by La Barbera et al. are at the 0.025 dex level, 
not yet accessible with our data. 

6. SUMMARY 

Using integral-field spectroscopy, we have looked at the aver¬ 
age stellar population gradients for a large sample of 100 early 
type massive galaxies. We are able to reach radii 2.5Re or ^ 
15 kpc. In keeping with previous results, we find no significant 
gradients in stellar population age nor [a/Fe] abundance ratios 
with radius, and gentle gradients in [Fe/H]. We thus confirm 
our previous result that the stellar populations in the outskirts 
of massive galaxies have sub-solar [Fe/H] but are enhanced in 
[a/Fe], suggesting th at the stars formed quickly and early, b ut 
in shallow potentials (iGreene et al.ll2013tlBender et alJl2015h . 

We examine the stellar population properties weighted to¬ 
wards 3-6 kpc, the ty pical sizes of massive g alaxy cores as ob¬ 
served at z Ri 2 (e.g.. Ivan der Wei et al.ll201^ . We see that age 
and [a/Fe] rise with increasing cr*, as we might expect if denser 
galaxies with higher cr* form earlier. We also find that at fixed 
physical radius, [Fe/H] and [C/Fe] correlate more strongly with 
Mfc than stellar velocity dispersion. In contrast, when look¬ 
ing at bins weighted towards ^ Re, we find no strong trends 
between abundances or abundance ratios and cr* or Mk- The 
average star as measured near the half-light radius in the most 
massive ellipticals is similar to the average star as measured at 
the half-light radius in galaxies of lower mass, as we might ex¬ 
pect if large galaxies grow via accreting smaller satellites. We 
do, however, still see a trend between stellar age and cr* even in 
the /?e-weighted bins. 

The gradients in [C/Fe] are similar to those seen in [Fe/H]. 
We suggest that the C comes mainly from mass loss in massive 
stars because th ere is not time to ge t it from intermediate-mass 
AGB stars (e.g. . [Graves et alJl2()07ll. Higher y ields due to mass 
loss from metal-rich stars (e.g.. lMaedeilll992h cause a pseudo¬ 
secondary dependence of carbon on Fe. In contrast, we see 
super-solar [N/Fe] that persists to large radius; the mismatch 
between [C/Fe] and [N/Fe] is a puzzle. [Ca/Fe] has solar values 
over the entire radial range that we observe (with the possible 
exception of a rise at large radius in the high dispersion bin), 
consistent with the idea that significant Ca is produced in Type 
la supernovae. 

Thanks to our relatively large sample, we are able to ex¬ 
amine trends in stellar mass at fixed cr*. At fixed cr*, we 
find marginal evidence that galaxies with higher stellar mass 
are more a-enhanced and [Ca/Fe] enhanced (because they are 
slightly [Fe/H] poor), suggesting a shorter timescale for star 
formation in more massive systems. 

Finally, we perform one of the most extensive spectroscopic 
studies of stellar population gradients as a function of group 
richness, while controlling for cr*. Overall, the differences in 
stellar population properties at large radius as a function of 
richness are very small, suggesting that internal properties like 
cr* deter mine stellar popu lation gradients (as in galaxy cen¬ 
ters; e.g.. IZhu et al.l 1201 ^. Galaxies in richer environments 
(/Vnei > 3) tend to be slightly older, slightly a-enhanced, and 
slightly [Fe/H] poor. We also see very slight trends towards 
shallower declines in [Fe/H] in lower-mass halos. When we re¬ 
strict attention to only central galaxies, these slight differences 
vanish, perhaps suggesting that they are driven by the fraction 
of satellite galaxies in the stacks. These trends are quite weak. 
Better statistics are needed to confirm them. Furthermore, it 
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Fig. 9. — Radial gradients in age, [Fe/H], [Mg/Fe], [C/Fe], [N/Fe], and [Ca/Fe] as above, but now in bins of group richness for so-called central galaxies (the 
most massive galaxy in the group, according the the Crook et al. 2007 catalog). As above, we divide the sample into those with three or fewer neighbors, and all the 
rest. We also control the distribution in cr* to match between the two bins. We fit the radial gradients with a power law of the form X = A log(R//? 3 ) + B for each 
stellar population parameter X, where is either 3-6 kpc or 1 - 1 .5Re . The fits to the rich (long-dashed lines) and poor (dotted lines) bins are shown here and in 
Tables 1 & 2. 
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will be quite interesting to c ombine our dynamical and stellar 
population information (e.g.. lRaskutti et ani2014t|Jimmv et^ 
l20lllNaabetai1l20l4 . 


By the end of the MASSIVE survey, we should roughly dou¬ 
ble the number of Mk < -25.3 mag galaxies in the sample, 
improving our ability to examine trends with environment and 
mass at fixed cr*. Our sample will be further complemented by 
ongo ing ambitious integra l-field gala xy surveys such as CAL- 
IFA (Isinchez et alJ|2012h . MaNGA dBundy et all 1201 5l) . and 


SAMI dCroom et ^ 20121) . When combined with our dynami¬ 


cal constraints on the total masses of these galaxies, these stellar 
population constraints on the age and the mass-to-light ratios of 
the galaxies will help address a number of pressing questions in 
galaxy eyolution, including the dependence of the initial mass 
function of cr*, the ratio of black hole mass to stellar mass at 
the high mass end, and the role of dark matter halo mass in the 
internal eyolution of massiye galaxies. 
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